On the half-metallicity of Co2FeSi Heusler alloy: an experimental and ab initio study 
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Co2FeSi, a Heusler alloy with the highest magnetic moment per unit cell and the highest Curie 
temperature, has largely been described theoretically as a half-metal. This conclusion, however, 
disagrees with Point Contact Andreev Reflection (PCAR) spectroscopy measurements, which give 
much lower values of spin polarization, P. Here, we present the spin polarization measurements 
of Co2FeSi by the PCAR technique, along with a thorough computational exploration, within the 
DFT and a GGA+U approach, of the Coulomb exchange U-parameters for Co and Fe atoms, taking 
into account spin-orbit coupling. We find that the orbital contribution (rrio) to the total magnetic 
moment (mr) is significant, since it is at least 3 times greater than the experimental uncertainty 
of niT- Account of ruo radically affects the acceptable values of U. Specifically, we find no values of 
U that would simultaneously satisfy the experimental values of the magnetic moment and result in 
the half-metallicity of Co2FeSi. On the other hand, the ranges of U that we report as acceptable 
are compatible with spin polarization measurements (ours and the ones found in the literature), 
which all are within approximately 40-60% range. Thus, based on reconciling experimental and 
computational results, we conclude that: a) spin-orbit coupling cannot be neglected in calculating 
Co2FeSi magnetic properties, and b) Co2FeSi Heusler alloy is not half- metallic. We believe that our 
approach can be applied to other Heusler alloys such as Co2FeAl. 

PACS numbers: 71.20.-b, 71.20.Be, 75.70.Cn 
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While Heusler compounds have been known for more 
than a hundred yearsi, they drew a remarkable amount 
of attention^Ti^ ever since the prediction by de Groot et 
ali^ in the early 1980's that some of these alloys would 
have a metallic band structure for the majority spin chan- 
nel and a semiconducting band structure for the minor- 
ity one, resulting in 100% spin-polarization (P) at the 
Fermi level. Such iiaif-metallic (HM) materials, with high 
values of P and Curie temperature (Tc), are excellent 
candidates for spintronic applications^iiiS (e-g-: magnetic 
random access memories (MRAM)ii utilizing the giant 
magneto-resistance spin- valve effecti^ in magnetic tunnel 
junctionai^ ). Specifically for Co2FeSi, high and low tem- 
perature magnetometry experimental^ showed that it is 
a Heusler compound with the highest magnetic moment 
((5.97 ± 0.05)/xb per unit cell, at 5 K) and the highest Tc 
(1100 K) among other Heusler alloys. 

Detailed computational studies have indicated that 
an orbital-dependent potential accounting for a mod- 
erate Coulomb-exchange interaction must be included 
in self-consistent calculations to simultaneously repli- 
cate both the experimental equilibrium lattice param- 
eter (5.64 A) and the measured magnetic moment of 
Co2FeSi alloy] ^^d^ These studies have also revealed that 
a total spin magnetic moment ~ 6 can be obtained 
only for the effective Coulomb-exchange interaction pa- 
rametersi^ [U^ff = U — J) falling within the ranges 
of 2.5-5 eV and 2.4-4.8 eV for the d-orbitals of Co 
and Fe atoms, respectively. Even though Co2FeSi alloy 
appears to be half-metallic only theoretically and only 



under stringent conditions on Ueff, it has been exten- 
sively referenced in the literature as suchjii"— This pre- 
diction, however, is at odds with several experimental 
measurements based on Point Contact Andreev Reflec- 
tion (PCAR) spectroscopj*i2i22ii^, which yield values of 
P -50%, far lower than 100%. 

The goal of this Rapid Communication is to recon- 
cile the results of computational predictions with experi- 
mental measurements for Co2FeSi. First, we present the 
results of our own PCAR measurements of P and com- 
pare them with those available in the literature. Second, 
we perform a thorough computational exploration of the 
Coulomb exchange ?7-parameter space for the 3d-orbitals 
of Co and Fe atoms in Co2FeSi, seeking the domain of 
parameters allowing to replicate the experimental mea- 
surements. 

The samples were prepared by arc melting of stoichio- 
metric amounts of the constituents in an argon atmo- 
sphere at 10 mbar. The polycrystalline ingots were 
then annealed in an evacuated quartz tube at 1273 K 
for 21 days. This procedure resulted in samples exhibit- 
ing the Heusler type L2i structure, which was verified 
by X-ray powder diffraction (XRD) using excitation by 
Mo Kqi radiation. Flat disks were cut from the ingots 
and polished for spectroscopic investigations of bulk sam- 
ples. X-ray photo emission (ESCA) was used to verify 
the composition and to check the cleanliness of the sam- 
ples. After removal of the native oxide from the polished 
surfaces by Ar+ ion bombardment, no impurities were 



detected with ESCA. 

For the PCAR measurements, niobium superconduct- 
ing tips are fabricated by electrochemical etching of 250 
fim thick niobium wire in a solution of HNO3, HF, and 
CH3COOH, with a mixing ratio of 5:4:1 by volume. The 
wire was kept at a positive potential with respect to the 
graphite counter electrode. The applied voltage was op- 
timized at ~ 8 — 10 V for the output current of ~ 800-0 
mA to get a sharp tip, see Fig. [1^) and b). Just before 
the measurements the tip was briefly dipped into the HF 
solution. A freshly etched superconducting Nb tip (bulk 
Tc ^ 9.3 K) was then mounted onto a shaft connected to 
a differential type screw that could be driven manually 
by 10 /im per revolution. For the low temperature mea- 
surements both the tip and the sample were immersed 
into a liquid He bath. The current- voltage (I — V) mea- 
surements were taken using a four-probe technique, with 
the differential conductance dl/dV obtained by standard 
ac lock-in detection at a frequency of 2 kHz within the 
temperature range 1.2 — 4.2 K. Typical results for nor- 
malized conductance G{V)/Gn as a function of voltage 
V are shown in Fig. [T];). At least 15 different junc- 
tions with the contact resistance 1 < Rc < lOOJl were 
measured and analized. To extract the values of spin po- 
larization, P, the conductance curves for each junction 
were fitted with the modified^! BTK model2^, using the 
value of Nb superconducting gap, A = 1.5 meV. As the 
values of P were found practically independent of the in- 
terface transparency Z, the values for individual contacts 
were averaged out. The red circles in Fig. [TJ;) represent 
the experimental data and the dashed black line is the 
fit. We find that the average value of spin polarization 
(P) =48±3%, Table U shows that our results are in good 
agreement with the other PCAR measurements available 
in the literature. 



TABLE I: Experimental spin polarizations of Co2FeSi by the 
PCAR technique. 
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PCAR on CoaFeSi 




Rcf. 


48 ±3 


bulk @1.2 K 




this work 


49 ±2 


thin films on MgO (001) <i 


M.2 K 


17 


57 ± 1 


bulk @4.2 K 




23 


59 ±2 


thin films on n-Ge(lll) <e 


4.2 K 


22 



In order to proceed, we need to calculate the spin po- 
larization, P„, that can be defined by2&"— : 



Pr.. 



X 100 



(1) 



where the averages are taken upon all the sheets of the 
Fermi surface, and the exponent n depends on the details 
of the experimental technique. For spin-resolved photoe- 
mission n = (Pq is the "static" or "intrinsic" spin polar- 
ization); n = 1 corresponds to experiments in the ballistic 
transport regime, whereas n — 2 describes experiments 
dominated by diffusive transport. Ideally the PCAR ex- 
periments arc done in the ballistic (Sharvin) regime, but 
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FIG. 1: (Color online) Superconducting Nb tip: a) shows a 
scanning electron micrograph of a tip (x300), and b) shows 
just the apex of a tip from a difi^erent angle and a greater am- 
plification (x 60000). c-d) Optimum normalized conductance 
curves for Nb/CoaFeSi contacts as a function of voltage with 
Nb superconducting gap, A = 1.5 meV and T= 1.2 K. The red 
circles represent the experimental data and the dashed black 
line is the fit. 



if the mean free path is smaller than the minimum size 
of the contact, they can only be performed in the dif- 
fusive regime<^i^ It is clear from Eq. ([T|) that while Pq 
can be directly calculated from the spin polarized den- 
sity of states (DOS), Pi and P2, in addition to the DOS, 
also require the respective Fermi velocities. Following 
the approach of Scheidemantel et al^, we calculated the 
Fermi velocity directly from the matrix elements of the 
momentum operator (instead of differentiating Bloch en- 
ergies j: with respect to k) according to: 
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(2) 



where i is the band index. These matrix elements can 
be readily generated using the optical package of the 
full potential linearized augmented plane waves (FP- 
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FIG. 2: (Color online) Spin polarizations (Pq, Pi and P2, see text for details), and spin (rria) and total (spin plus orbital: 
rriT ~ rus + ruo) magnetic moment for selected values of Uco (from left to right: 0, 0.4, 0.8, 1.6, 2.8, 3.6 eV) and Upe in the 
range 0-6 eV. The gray shaded areas mark the experimental ranges for total magnetic momenli^^: mr = 5.97 ± 0.05 and the 
spin polarization (see Table[I]): 45 < P < 61. 



LAPW)^"^ WIEN2k code^, and further used as an in- 
put in Eq. ([T|) to calculate Pi and P2- We tested this 
procedure for pure bcc Fe and fee Ni. As we can see 
from Table our results are in satisfactory agreement 
with other calculations and experimental data available 
in the literature. 



TABLE II: Our calculated spin polarizations (in %) for bcc Fe 
and fee Ni, against others' calculations (figs. 2-3 in "Refi^, 
"Ref.^) and experiments ("Ref.^^, ^Ref.^^). 







bcc Fe 




fee Ni 


Ours 


Others 


Exp. 


Ours 


Others Exp. 


Po 


51 


58°, 59" 




-80 


-81° 


Pi 


39 


39°, 33" 




-45 


-48° 


P2 


37 


33°, 21" 




5 


2° 


p 






44±3\ 40-48^= 




40-47.5^^ 



We carried out all our calculations with the WIEN2K 
code, using the generalized gradient approximation 
(GGA) in the formal parameterization scheme of Perdew- 
Burke-Ernzerhof (PBE)4- and the experimental lattice 
parameter of 5.64 A for the cubic Co2FeSi (crystallo- 
graphic details of its structure can be found in Ref4i). 
Muffin tin radn (Rmt) of 2.32, 2.32, and 2.18 atomic 
units were used for the Co, Fe, and Si atoms, respectively. 
The R]^,Jx ' ^max product, where R]^,Jx is the smallest of 
all muffin tin radii, and K^ax is the plane wave cut-off, 
was set equal to 7 (implying a plane wave expansion cut- 
off of ^ 9.10 Ry), and the energy threshold between core 
and valence states used was of —6.0 Ry. Integration in 
the irreducible Brillouin zone was carried out over 641 
k-points and convergence was set to simultaneously be 
better than 10"'' (Ry) for the total energy and 10^^ (au) 
for the total electronic charge. Spin-orbit coupling was 
included in all calculations in order to obtain not only 
the total spin, but also the total orbital moment. Within 
the LDA(GGA)-hU schemei^ we have modified the 3d 
orbitals of Co and Fe atoms, and explored a square mesh 
of 16x16 points in the Uco-U pe space (resulting in a total 



of 256 self-consistent calculations), with < Uco < 6 eV 
and < [/fe < 6 eV. Since it has been shown^ that 
the choice of different "flavors" of the double counting 
correction of the -|-U method can be critical, we tested 
two of them: one, the so-called "SIC" (self-interaction 
correction), introduced by Anisimov et al.— and the 
"AMF" (around mean field), introduced by Czyzyk and 
Sawatzk^. The latter yielded the total magnetic mo- 
ments (spin plus orbital terms) lower than the lowest 
bound of the experiments (5.92 hb)^ up to L/ ~ 5 eV, 
so we found it inadequate for modeling Co2FeSi. Thus, 
all the results reported here were obtained using the SIC 
flavor of the method. 

Representative sampling of our results is shown in Fig. 
[51 where, the total spin moment (tos), the total magnetic 
moment (tot, i-C-, total spin (m^) plus total orbital (toq) 
moment), and the spin polarizations Pq, Pi, P2 are pre- 
sented as a function of Uco and Upe- It is clear that for 
a system where spin-orbit coupling is important, such as 
the case for Co2FeSi, the spin-orbit interaction induces a 
strong orbital component in the tot, which turns out to 
be significant and cannot be neglected: it shifts to lower 
energies, and simultaneously narrows down, the range of 
U Fe that yields a total magnetization within the experi- 
mental margin of error, tot = 5.97 ±0.05^ (see the gray 
shaded area in the bottom panels of Fig. [5]). One can 
also see from Fig. [2] that for Uco > 3.6 eV none of the 
C/Fe-values would yield the value of mx within the ex- 
perimental range. The total spin moment is also shown 
in all of the bottom panels (it is fixed at 6 /is beyond 
a certain value of Upe)- The analysis of all our calcula- 
tions allows us to produce Figs. |3h.) and b), that show 
the values of the U-parameter (areas within the closed 
loops in the figure) for which the calculated total mag- 
netic moment and spin polarizations {Pq, Pi, or P2) are 
within the experimental range. Finally, Fig. [3Ji) shows 
the values of U for which our calculations simultaneously 
agree with the results of both: magnetomctry and PCAR 
spectroscopy (the latter is taken over the entire range, 
45% < P < 61%, defined by our measurements and oth- 
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ers, see Table |T]). Except for the lower area {U ~ 0.4 eV , 
which corresponds to P2)j all other values correspond 
to Pi, i.e., the expected ballistic transport between the 
superconducting tip and the sample in the PCAR ex- 
periments with Co2FeSi. Thus, in view of our findings, 
the results of magnetometry measurements indicate that 
Co2FeSi is actually not a half-metal, and the determina- 
tion of P using PCAR spectroscopy is fully compatible 
with this prediction. 
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FIG. 3: (Color online) The area within the closed loops cor- 
respond to values of U for which a) the calculated total mag- 
netic moment, mr, agrees with experiment; b) the spin po- 
larizations (Po, Pi, or P2) are within the experimental PCAR 
values, c) The four loops of a) and b) are superimposed (omit- 
ting the symbols) to aid the eye. d) Intersection of the areas 
of a) and b), i.e., the values of U for which both calculated 
rriT and P values fall within experimental values. 

Chalsani et al^ mentioned that the results of the 
PCAR technique can be affected by a number of fac- 
tors, such as the geometry of the contact and interac- 
tions between the sample and the tip through surface 
states; however, ab initio calculations by Khosravizadch 



ct al4^ showed that the surface of Co2FeSi is also non 
half-metallic, and the loss half-metallicity of Co2FeSi, 
which we have found, cannot be attributed to the sur- 
face effects. Moreover, our results for a perfect bulk ma- 
terial would persist, regardless of whether or not anti-site 
defects, finite tcmperaturo^i, and crystallographic disor- 
der——, could induce additional loss of half-metallicity 
by the appearance of spin states in the minority spin 
channel. The orbital moments we calculated in Co2FeSi 
depend on the U-parameter, ranging from 0.145 /xb for 
Uco = Upe = 0, to 0.482 HB for Uco = Upe = & eV. One 
key idea presented here is that even the lowest orbital mo- 
ment obtained (toq = 0.145 mus) is almost three times 
greater than the experimental uncertainty (0.05 /is) in 
the magnetometry measurements of the total magnetic 
moment per formula unit^^ {mx ~ 5.97 ± 0.05/iB), as 
Fig. [2] clearly demonstrates. In addition, the total spin 
moment as is shown in Fig. [5] (in all the bottom panels) 
is fixed at 6 /.is beyond a certain value of Upe- It is also 
seen that for Uco ~ 3.6 eV there is no Vpe value that 
would yield an within the experimental range. Our 
calculations of mo, are also supported by the systematic 
study by Galanakis^. performed on nine fuU-Hcuslers 
alloys, and while the author does not study Co2FcSi, he 
presents the results on similar compounds, Co2FeAl and 
Co2MnSn, citing for the former the greatest total orbital 
moment calculated in that work to be 0.149 /.is, which is 
in excellent agreement with our 0.145 \ib value. 

In summary, our calculations provide a strong evidence 
that the orbital component of the total magnetic mo- 
ment in Co2FcSi cannot be neglected. By taking it into 
account we identify the ranges of the JJ-parameters com- 
patible with both the magnetometry and PCAR (ours 
and others) measurements. Based on the range of the 
U-parameters, we conclude that Co2FeSi is not a half- 
metal. We believe that our approach will be applicable 
to other compounds similar to Co2FeSi, such as Co2FeAl 
and Co2MnSn. 
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